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, The magnetic field strengths of most millisecond pulsars(MSP) are 

about 10 8-9 Gauss. The accretion induced magnetic field evolution sce- 
5J" 1 , nario here concludes that the field decay is invesely related to the accreted 

C/3 ' mass and the minimum field or bottom field stops at about 10 8 Gauss if 

— |- , accreted with the Eddington accretion rate, which is proportionally re- 

fN| ' lated with the accretion rate as M 1 / 2 . The possibility of the low field 

~ 10 7 Gauss MSPs has been proposed for the future radio observation. 



As known, the low magnetic field (MF) and fast spinning of MSPs are 
ascribed to the mass accretion in the binary system. The accretion induced 
neutron star magnetic field (NSMF) decay models have been proposed and paid 
much attention by a couple of authors (Romani 1990; Urpin & Geppert 1995; 
Cheng & Zhang 1998; Brown & Bildsten 1998; Payne & Melatos 2003), which is 
based on the observational evidence that NSMF decays in the binary accretion 
phase (Taam & van den Heuvel 1986; Shibazaki et al. 1989). Moreover, Van 
Q_i' den Heuvel and Bitzaraki (1995) discovered that the increased amount of mass 

accreted leads to decay of the NSMF, and the 'bottom' field strength of about 
| 10 8 Gauss is also implied. 

Based on the accretion induced field decay in the slab geometry approxima- 
tion proposed by Cheng & Zhang (1998), we extend it to the the spherical geome- 
try structure of the accreted neutron star illustrated in Fig.l of paper by Cheng & 
Zhang (1998). We assume that the magnetic field lines are frozened in the entire 
NS crust with homogenous average mass density. Under the condition of incom- 
pressible fluid approximation as well as the constant crustal volume assumption, 
the piled accreted materials will arise the expansion of the volume of the polar 

zone by SV P = — A p 5H, where A p is the surface area of the polar zone, p is 

the average density of the crust, M is the accretion rate, 5t a particular accretion 
duration, H the thickness of the crust and 5H/H the fraction of the thickness 

dissolving into the core which is defined by ^jf = where M cr = 4nR 2 pH is 
the crust mass. The expansion of the polar zone will dilute the magnetic flux 
density if the conservation of the magnetic flux is preserved, so, 5(BA p ) = 0, 
and the area A p of the accretion polar patch can be accurately expressed as 
follows(Shapiro & Teukolsky 1983), A p = 2irR 2 (1 - cosO c ) , sin 2 c = where 
C is the angle between field line of the star surface and polar axis, Ra is the 
Alfven radius, R A = 3.2 x 10 8 (cm) M~ 2/7 p% 7 (^L)- 1 / 7 , where M is the neu- 
tron star mass, Myj is the accretion rate in units of W 17 g/s, and p^o is the 
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magnetic moment in unit of 10 30 Gem 3 . Therefore, using the specified relation 
5V P = H5A p and the equation of conservation of magnetic flux, we get the fol- 
lowing relation, 6V P = —V p ^-. Connecting the equations above, we obtain MF 
evolutionary equation as, A P 5B /[(2itR 2 — A p )B\ = M5t/M cr . Considering the 
initial condition B(t=0)=So, we have, 

B/B f = {1 - [Cexp(-^g) - 1] 2 }" 7/4 , C = 1 + v/l - 

AM = Mt and Bf is the magnetic field defined by the Alfven radius matching 
the star radius, i.e., Ra{Bj) = R, which is also the minimum field strength, or 
named as the bottom field (van den Heuvel & Bitzaraki 1995), 

B f = 1.32 X 108(^)1/2(^)1/4^5/4 (G)j 

where Mecl is the Eddington accretion rate and Rq is NS radius in units of 
10 6 cm. If AM = 1CT 5 ~ 1(T 4 M , MF evolutionary equation can be sim- 
plified as the following approximated form, B = Bq/(1 + 7^), where tub = 

±(B f /B ) A/7 M cr ~ 5 x lO~ 4 (M )(M cr /O.2M Q ). This is just the same form as 
the empirical formula of accretion induced field decay proposed by Shibazaki et 
al. (1989). 

The main conclusions are summaried in the following. NSMF decays in the 
binary accretion phase, which is inversely correlated with the mass acreted from 
the companion. The bottom field of MSP is determined by the condition of that 
the magnetosphere radius equals the star radius, which is proportionally related 
to the accretion rate as M 1 / 2 (White & Zhang 1997). The bottom field of Z- 
source (Eddington accretion rate) is about ~ 10 8 Gauss, and the bottom field 
of Atoll-source (1% Eddington accretion rate) will correspond to ~ 10 7 Gauss. 
The so low field ~ 10 Gauss MSPs have not yet been discovered from the radio 
observation. The final state of NS is constrained by the system parameters, 
such as star radius, star mass and mass accretion rate, which is nothing to do 
with the initial field and initial period, so this is why the almost homogeneous 
field distributions of MSPs from ~ 10 8 (G) to ~ 10 9 (G) do not follow the field 
distributions of the normal PSRs from ~ 10 11 (G) to ~ 10 15 (G). 
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